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ABSTRACT 

Context. Based on the far infrared excess the Herbig class of stars is divided into a group with flaring circumstellar disks (group I) 
and a group with flat circumstellar disks (group II). Dust sedimentation is generally proposed as an evolution mechanism to transform 
flaring disks into flat disks. Theory predicts that during this process the disks preserve their gas content, however observations of 
group II Herbig Ae stars demonstrate a lack of gas. 

Aims. We map the spatial distribution of the gas and dust around the group II Herbig Ae star HD 95881. 

Methods. We analyze optical photometry, Q-band imaging, infrared spectroscopy, and K and N-band interferometric spectroscopy. 
We use a Monte Carlo radiative transfer code to create a model for the density and temperature structure which quite accurately 
reproduces all the observables. 

Results. We derive a consistent picture in which the disk consists of a thick puffed up inner rim and an outer region which has a 
flaring gas surface and is relatively void of 'visible' dust grains. 

Conclusions. HD 95881 is in a transition phase from a gas rich flaring disk to a gas poor self-shadowed disk. 
Key words, star formation, protoplanetary disks, Herbig Ae stars, HD 95881 



[ 1. Introduction 

Herbig Ae (HAe) stars are known to have gas-rich, dusty disks 
that are the remnant of the star formation process. These disks 
are most likely the sites of ongoing planet formation. The pro- 
cesses leading to and associated with planet formation modify 
both the composition and the geometry of the disk. Grain growth 
and grain settling are expected to result in large spatial variation 
of the grain size distribution and the gas to dust mass ratio within 
the disk. The gravitational interaction of proto-planets with the 
i disk can create gaps/holes. Also the gas and the dust chemistry 
' is expected to vary spatially. In order to understand planet for- 
mation, it is thus important to establish the spatial distribution of 
gas and dust in protoplanetary disks independently. 
[ Observationally, the Spectral Energy Distributions (SEDs) of 
1 HAe stars have been divided into two groups, that reflect dif- 
ferences in the slop e of the mid-IR (10-60 jum) spectral range 
(iMeeus et al 1 l20()lh . Group! sources have red SEDs, while 
group II sources have blue SEDs. These diff'erences can be in- 
terpreted in terms of the geometry of the disk. The direct irra- 
diation of the inner rim of a disk with an inner hole causes it to 
be puffed up (Dullemond et al., 2001). This puffed up inner rim 
casts a shadow, and only the outer disk surface regions emerge 



from the shadow and receive direct stellar light. Depending 
on the dust opacity, some disks may never eme rge from the 
shadow of the inner rim dDuUemond & Dominikl l2004i) . This 
provides an elegant explanation for the observed two types of 
SEDs: group I sources being flaring, and group II sources self- 
shadowed. This interpretation has been confirmed using spa- 
tially resolved mid-infrared (IR) im aging with the Ver y Large 
Telescope Interferometer (VLTI; e.g. lLeinert et al.l&004l) . 

A dilference between group I and group II sources was also 
found for the strength of the mid-IR emissio n bands attributed to 
Polycyclic Aromatic Hydroc arbons (PAHs; iMeeus et~al ] l200Tt 
lAcke & van den Ancken2004l) : flaring disks tend to show strong 
PAH emission while self-shadow ed sources show weaker or no 
PAH emission (lAcke et al.l2010h . A similar d ilference was foun d 
for the strength of flie [OI] 6300 A line (lAcke et aI.L l2005h . 
However, there is significant scatter in these trends (see below). 
Both the PAHs and the [O I] line strength probe the gas in the up- 
per disk layers, and both require direct irradiation of the disk sur- 
face by stellar photons to be excited. PA Hs mainly pro be the disk 



on scales of several tens to 100 AU (e . g.lvan Boekel et al ^ool 
iLagage et afll2006t lOeers et 
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I.e. 

similar scales as the dust continuum emission in the 10-60 pm 
wavelength range. These observations suggest that the spatial 
distribution of gas in group II sources is different from that of 
group I sources: apparently, in group II sources the gas in the 
surface of the outer disk does not receive direct stellar photons. 

In a theoretical study, iDuIIemond et al.l (l2007b show that 
for disks in which the dust settles but the scale-height of 
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Table 1. Characterizing parameters of HD 95881. The sec- 
ond C ol, lists parameters as given by lAcke & van den Anckerl 
( |2004|) . where SFR stands for the star formation region. New 
values for some of the parameters are given in the third Col. 



Parameter 


Value 


New 


Right Ascension 


ll*" 01'" 57^62 




Declination 


-72° 30' 48'.'4 




Spectral Type 


A2III/IVe 




T-eff [K] 


8990 




Distance [pc] 


118+12 


170+30 


Luminosity [Lq] 


6.9±1.0 


15.4±6 


Am [mag] 


0.25 




Radius [Rq] 


1.1±0.1 


1.6+0.3 


Mass [Mq] 




2.0+0.3 


SFR 


Sco OB2-4? 


? 


Group 


Ila 





we consider the determination of T^g as reliable, we adopted 
a luminosity of 15.4 assuming the ZAMS luminosity from 
Mevnet et al. ( 1994) at the given T^g. The corresponding ZAMS 
mass was adopted as the stellar mass. From the new luminosity, 
new estimates for the radius and distance follow. Table [T] lists 
the basic characteristics and our new estimates. To double check 
the new distance we looked at the Tycho parallax, which gave 
a lower limit of 80 pc. We also consulted the Hipparcos cata- 
log for spectral types and (B - V) photometry of stars within 1° 
from HD 95881 and com pared the distaii ce to the color excess 
E(B - V) for the region dde Zeeuw et"al]|1999l) . For HD 95881 
the color excess of E{B - V) - 0.08 leads to an upper limit 
on the distance of ~200 pc. The derived distance interval of 80- 
200 pc is consistent with both distance estimates, but because of 
the argument given above we adopt a distance estimate as de- 
rived from the ZAMS of 170+30pc. 



the gas does not change, both group I and group II sources 
should show prominent PAH emission, contrary to the ob- 
served trend. However, some disks classified as group II 
sources (i.e. with a self-shadowed dust geometry) are ob- 
served to show prominent PAH emission an d [OI] line emis- 
sion; examples are HD 98922 a i id HD 95881 dAcke et all l2005t 
lAcke & van den Anckeiil2004 . IVerhoe ff et al.l (l2010h show that 
for HD 95881 the PAH emission is extended at a scale sim- 
ilar to those of group I sources. iFedele etakl (l2008l) and 
Ivan der Plas et al.l ( 20081) studied three HAeBe stars and found 
that the gas and dust in their disks has a different spatial dis- 
tribution. In the case of the group II source HD 101412, PAH 
and [O I] emission were detected and found to be more extended 
than the dust continuum at 10 /im. These observations suggest 
that disks exist in which the dust has settled but the scale-height 
of the gas is still large enough at several tens of AU distance 
from the star to produce substantial PAH and [01] line emis- 
sion. Such disks may provide important clues as to how gas-rich 
disks evolve to gas-poor debris disks. 

In this paper, we study the spatial distribution of the gas and 
dust in the disk of HD 95881. This star was part of a larger 
study of spatially resolved mid-IR spectroscopy of HAe stars 
dVerhoeff et al.. 2010) . We use optical spectroscopy of the [OI] 
line, the SED, infrared spectra as well as near-IR and mid-IR in- 
terferometric observations to constrain the geometry of the gas 
and dust in the disk. We use a hydrostatic equilibrium disk model 
to fit the SED and compare the predicted spatial distribution of 
the near-IR and mid-IR emission of the best fitting disk model 
to the interferometric observations. We find convincing evidence 
that the dust in the disk of HD 95881 has setded but that the gas 
still has a significant scale-height. We derive an estimate on the 
total disk mass by fitting the strength of the PAH bands. 



2. Stellar parameters 

There is little information available on HD 95881 in the litera- 
ture. From optica l photometry the spectr al type was determined 
to be A2III/IVe (Houk & Cowlevlll975|)^ which was translated 
into a n effective temperature of 8990 K dXcke & van den Ancker 
12004 . The distance was established by a relatively uncer- 
tain association with the star formation region Sco OB2-4 
dAcke & van den Anckeiil200l. The distance to this region was 
previously determined bv'd e Zeeuw et al.l (119991 When we take 
these parameters to pinpoint HD 95881 in the Hertzsprung- 
Russell-diagram, we find that it is situated to the left of the 
Zero Age Main Sequence (ZAMS), which is unphysical. Since 



3. Observations 

3.1. AMBER 

Spectrally dispersed K-band interferometric observations of 
HD 95881 were obtained with VLTI/AMBER on the UT1-UT3- 
UT4 baseline setting during the night of February 4, 2007. The 
weather conditions were excellent, with the optical seeing as 
low as 0.5". The data have been reduced according to t he stan- 
dard AMBER data reduction procedure described in Tatull i etaP 
(2007). AMBER observations of a standard star (HD 100901, 
KO/KIIII) were performed directly before the science measure- 
ment and used to calibrate the raw visibilities. 



3.2. VISIR imaging 

Q-band imaging data were obtained on the 14th of December 
2005 using t he VLT Imager and Spectrometer for mid-Infrared 
(VISIR; see iLagage et al. |2004|) . Standard "chopping and nod- 
ding" mid-infrared observational technique was used to suppress 
the background dominating at these wavelength. The Q2 filter 
was chosen, which has a central wavelength of 18.72yum and a 
half-band width of 0.88 /im. The pixel field of view was 0.075" 
and the orientation was standard (north up, east to the left). The 
total integration time spend on HD 95881 was 680 s. The ref- 
erence star, HD 1024601 (10.6 Jy in the Q2 filter), was chosen 
from the database of VISIR standard stars based on criteria of 
flux and distance on the sky to the target. 456 s of integration 
time was spent on HD 102461. The airmass of both sources was 
always below 1.5. The sensitivity achieved was estimated to be 
57 mJy / 10 cr/lh;. The optical seeing was moderately fluctuating 
in the range 0.75-0.85", corresponding to an estimated Q-band 
seeing of <0.4". 



3.3. VISIR spectroscopy 

Long slit N-band spectra were obtained with VISIR in the 
low resolution (LR) mode. A sample of 17 HAe stars was ob- 
served u nder the VISIR GTO program on circumstellar disks 
(IVerhoeffet al.ll2010l) . During the nights of December 16, 2005 
and April 8, 2006 we observed HD 95881. Standard chopping 
and nodding along the slit was used with a chopper throw of 8", 
a slit-width of 0.75", and a standard North-South orientation. 
After and before the science measurements bright standard stars 
were observed. The airmass of the observations was ~1.6 and the 
optical seeing was around 1 .0" during the first night and around 
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1 .2" during the second night, corresponding to an estimated N- 
band seeing of <0.65". 

3.4. MIDI 

HD 95881 was observ ed with the MID-in fmred Interferometric 
instrument (MIDI; see Leinert et al l200l . the lOyum instrument 
of the VLTI, during the night starting June 6, 2004, as part of the 
"science demonstration" program. The UT1-UT3 baseline was 
used, resulting in a projected baseline length of 65 m at a position 
angle of 81° east of north. The grism was used to spectrally dis- 
perse the signal at a resolution of R = Ajt^A - 230. This ensures 
that relatively narrow "dusty" emission features such as those of 
crystalline silicates and PAHs are well resolved spectrally. The 
seeing was constant at a value of 0.6" and the atmospheric trans- 
parency was excellent. 

We performed observations in "High-Sens" mode, i.e. we 
took an interferometric measurement combining the light from 
both telescopes, and subsequent photometric measurements of 
the signal from one telescope at a time. Stars of known bright- 
ness and angular diameter were observed, using the same pro- 
cedure, for photometric calibration and monitoring system co- 
herence losses (interferometric "transfer function"). We found 
our interferometric measurement of HD 95881 to be of signifi- 
cantly higher quality than the corresponding photometric mea- 
surements, and chose to directly calibrate the correlated flux 
rather than converting to interferometric visibility by division 
over the photometry. 

3.5. Additional data 

A low resolution Spitzer-IRS spectrum is used to com- 
pare and flux-calibrate the VISIR spectrum. For the ac- 
quisition and red uction of the Spitzer data we refer to 
Juhasz et"al] ( 20101). Photomet ric data points were taken from 



AMBER 



1.0 



Acke & van den Anckerl (l2004l) . Together with the Spitzer spec- 



trum they make up the SED, which allows us to constrain the 
disk geometry (see Sect.[5l). Finally, [ O I] data were taken from 
an earlier study, see lAcke et al.l (l2005h . 



4. Analysis 

4.1. AMBER 

A good fit (;^-^/v = 1.8) to the AMBER data could be obtained 
with an inclined uniformly emitting ring surrounding a point- 
like source representing the star (Fig.[T]i. The inner and outer 
radius of this ring are 2.2+0.2 mas (0.37 AU) and 2.7+0.3 mas 
(0.46 AU). The disk inclination is 60°+10° and its position an- 
gle is 102° +2°. The mean radius of the ring is 2.4 mas (-0.4 AU). 
This is similar to the values found by 'Ei sner et"al] ( |2004|) for a 
group of Herbig Ae/Be stars. The fraction of the K-band flux 
that is ascribed to the star by the fit of a Kurucz model to the 
optical photometry (~20%) is not exactly equal to the AMBER 
estimate of 32+4%. It may be that this discrepancy is due to 
the presence of an unresolved disk component, e.g. hot gas well 
within the dust sublimation radius, that contributes to the K-band 
emission. Gas emission from a very compact region around the 
star, where no dust species can survive, ha s been claimed for 
quite a nurnber o f pre-main-sequen c e stars (Eisner et al. 200 
iKraus et al.' 2008; Ack e et alJ l2008t iTannirkulamTt all l2008 
and Jsella et al. 2008). However, refractory dust has also been 
sugge sted to be present in the innermost regions dBenistv et alj 
I2009h . 



0.4 



0.2 



(36m, 65") 



_ (85m, -58°) 
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2.1 2.2 
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Fig. 1. The AMBER K-band visibiHties of HD 95881. The pro- 
jected length and position angle of the baselines is indicated 
on the left. Overplotted is the best-fit ring+point-source model 
(sohd fines; see text). 



1.10 




-20 
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Fig. 2. The continuum-norma Uzed [O I] 6300 A emission line of 
HD 95881 dAckeet al.ll2005h . The velocity axis has been cen- 
tered on the centroid position of the feature, i.e. the radial veloc- 
ity of the central star. Line profiles corresponding to a surface 
brightness proportional to r^^ are overplotted. Red represents a 
model with an outer radius beyond 3 1 AU; yellow corresponds 
to the model with an outer radius equal to 1 1 AU. See text for 
details. 



4.2. [01] data 



lAcke et aP (l2005h have investigated the [O I] 6300 A emission 
line in a large sample of Herbig stars. The authors argue that 
the emission is non-thermal and originates from the disk sur- 
face of a flared circumstellar disk. This is consistent with the 
large detection rate of [O I] emitters among the group I sources 
in the sample. Roughly half of the group II sources, however, dis- 
play the [OI] 6300 A line in emission as well, albeit less strong. 
HD 95881 is one of these targets. 

In Fig.|2l the [O I] 6300 A line profile is shown. We have fit- 
ted a generic model to the data, assuming that the intensity drops 
off with radius as a power law and that the disk is in Keplerian 
rotation. We adopt the disk inclination derived from the AMBER 
fit. An acceptable fit (reduced = 3) was achieved with a power 
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Fig. 3. The SpitzerLR spectrum of HD 95881 (SNR ^ 500). The 
left rake demonstrates the presence of the PAH bands at 6.2, 
7.9, 8.6, 11.2, 12.7, and 16.4 //m. The right rake points at the 
forsterite bands at 11.3, 16.2, 19.5, 23.5, and 33.5 /zm. The red 
square is the photometry from the VISIR Q-band imaging (see 
Sect.EHl. 



HD95881 Q2 filter 




0.8 1.0 1.2 
Distance to the star [arcsec] 



1.6 



Fig. 4. Surface brightness levels as a function of distance from 
the star in the Q2 image after point source subtraction. The error 
bars give the 3 cr uncertainty. 



index of -3.0+0.2, an inner radius of 0.9±0.2 AU and an outer 
radius larger than 11 AU. The best fit (reduced = 1.5) has 
an outer radius beyond 3 1 AU. In the figure we show the best- 
fit profile, as well as the profile that corresponds to the model 
with an 1 1 AU outer radius. Decreasing the outer radius of the 
model to even lower values would further deteriorate the fit qual- 
ity. We can therefore conclude that a significant fraction of the 
[O I] emission comes from a region at tens of AU from the star 
Note that our model does not account for the 15% of [O I] flux 
which is emitted in the blue- and red wings of the profile. In these 
regions, close to the star, the power law approximation fails. 

4.3. Spitzer 

The Spitzer-IRS spectrum is given in Fig.|3] It displays a rich 
mineralogy, there is emission of amorphous and crystall ine sili- 
cates as well as emission of various PAH bands. See Juha sz et al.l 
(1201 0) for a detailed discussion of the mineralogy. The blue 
slope at longer wavelengths is typical for group II sources. 

Compared to other HAe stars HD 95881 appears to be a 
very ordinary PAH emitter We have verified this by comparing 
PAH band-strength ratios. We have ta ken flie 6.2, 7.7 , 8.6, 11.2 
and 12.7 A/m PAH band strengths from lVerhoeff et afl (120 10) and 
lAcke et aTl ( 12010) and looked at the ratios of their continuum 
subtracted and integrated strengths. The only aspect that makes 
HD 95881 stand out sUghtly is the 8.6yum PAH feature, which is 
relatively strong. 

4.4. VISIR imaging 

A dedicated data reduction pipeline was used for the imaging. 
It features a comprehensive set of methods to correct for in- 
strumental signatures such as det ector striping or backg round 
low-frequencies excess of noise (iPantin et al.1 120081 12009). A 
photometric analysis gives an integrated flux for HD 95881 of 
6.0+0.5 Jy in the Q2 filter (18.72jum). This is in fair agreement 
with the Spitzer data (see Fig.O. Since the emission of the star 
is negligible at this wavelength, this flux can be attributed to 
the disk. We searched for an extended emission component us- 
ing PSF subtraction at a sub-pixel (1/10) precision level. The 
PSF was derived from the observation of the standard star The 
resulting residuals display an excess of signal that is roughly 



axi-symmetric and decreases gradually as a function of distance 
from the star. 

Assuming that the disk has the geometric parameters given 
in table|2] we divided the residuals in a series of concentric el- 
lipses having a separation of 0.075" along the semi-major axis 
of the disk. We assumed an aspect ratio corresponding to a flat 
disk inclined at 55°, and a position angle of 103° from North. 
We used these elliptic annuli to numerically estimate the surface 
brightness distribution and the corresponding uncertainty levels. 

Our statistical analysis confirms we have a significant detec- 
tion since for distances in the range 0.4-0.85" from the star the 
average values in each of the elliptical annuli are well above the 
detection limits set at 99% confidence. The measured average 
surface brightness in elliptic annuli as a function of the semi- 
major axis distance to the star are plotted in Fig.|4] The 3 cr un- 
certainties are displayed as bars on the plot. The total flux in the 
resolved component is 0.7+0.1 Jy. 

4.5. VISIR spectroscopy 

The general reduction and analysis strateg y of the sp ectroscopic 
VISIR data is described in detail in Verhoeff et al.( (i2010). For 
HD 95881 the telluric correction was done by means of an air- 
mass interpolation of two calibrators. The observation of April 8, 
2006, appeared to be suffering strong atmospheric residuals. The 
spectrum of December 16 2005 was thus chosen as most rep- 
resentative. The missing observation of the 9.8 fim setting was 
replaced with the poor one from April 8, 2006. We scaled the 
VISIR spectrum to the Spitzer spectrum at 10.6/im using an 
average factor for the different settings of 1.16. The resulting 
spectrum in Fig. |5] has a SNR of ~300. The agreement with the 
Spitzer spectrum is encouraging. The slight deviation observable 
just left of the ozone band at 9.6 //m is typical for the quaUty of 
the data taken on April 8, 2006. 

The FuU Width at Half Maximum (FWHM) of the spatial 
emission profile of the target was determined by performing a 
Gauss-fit in 32 merged wavelength bins. Comparison of the sci- 
ence signals with the PSF shows that the target is unresolved in 
the continuum. After quadratic subtraction of the PSF and av- 
eraging over the median values of all measurements we find a 
three sigma upper-limit to the FWHM extent of the continuum 
emission region of < 0.46", which corresponds to < 79 AU at the 
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Fig. 5. The top panel shows the VISIR N-band spectrum of 
HD 95881 (black). Overplotted are the Spitzer spectrum (thick 
green), the correlated flux as measured with MIDI (dotted blue), 
and the difference between the Spitzer and MIDI spectra (or- 
ange). The Spitzer spectrum represents the entire disk, the MIDI 
correlated flux spectrum represents the inner disk (r < 2 AU), 
and the diff'erence spectrum represents the outer disk (2 < 
r < 100 AU). The bottom panel shows the normalized contin- 
uum subtracted silicate feature (same color coding). The similar 
shape of the inner and outer disk silicate feature indicates a sim- 
ilar composition and grain size distribution. 



adopted distance of 170 pc. However the science signal displays 
a relative increase at 8.6 and 11.2//m and an upturn to the left 
of ~8 jum, which are exactly the wavelengths at which the PAH 
molecules have emission features (see Verhoeff et al. 201^. We 
checked the significance of these FWHM features with respect 
to pixel-to-pixel variations and concluded that the PAH emission 
is significantly more extended than the continuum. 

In order to estimate the spatial extent of the PAH emission 
we measured the spatial emission profile at the peak wavelengths 
of the PAH bands and we subtracted the spatial emission profile 
of the continuum contribution. This continuum profile was de- 
termined by interpolating the intensities and spatial profiles ad- 
jacent to the PAH bands. The resulting observed spatial profile 
of the PAH emission was Gaussian fitted to obtain the FWHM. 
Finally, the instrumental width (i.e. the PSF) was quadratically 
subtracted to obtain a measure for the intrinsic extent of the PAH 
emission. We found FWHM values of 0.34"!j;:j;^ and 0.39"!JS:[|^ 
for the 8.6 and 1 1 .2 jum PAH bands respectively, which results in 
absolute sizes of 58 and 66 AU. In a Gaussian distribution of the 
PAH emission this would mean that 99% is confined in a radius 
of ~100 AU. Note that this is a conservative estimate of the PAH 
emission scale since the PAH surface brightness is expected to 
fall off with distance from the star together with the stellar flux 
as 1/r^. 



4.6. MIDI 

In Fig.|5]we compare the spectrum in correlated flux as seen by 
MIDI, i.e. the total flux times the visibility, to the spectra ob- 
served by VISIR and Spitzer. The correlated flux spectrum is 
dominated by the central few AU of the disk, the Spitzer and 
VISIR spectra probe the entire disk. Note the difference in the 
strength of the 8.6 and 1 1 .2 /vm PAH bands in the spectra. These 
bands are prominent in the total flux spectra, but essentially ab- 
sent in the correlated flux spectrum. The PAH emission region 
is apparently outside of the disk region probed by MIDI. This 
shows by direct measurement that the PAH features arise at 
scales much larger than ~2AU in the disk ofHD 95881 . To stress 
this point we also plotted the difference between the Spitzer and 
the MIDI correlated flux spectrum in the bottom panel of Fig.|5] 
This difference spectrum is dominated by the emission of the 
outer disk (2 < r < 100 AU) and shows very distinct PAH fea- 
tures. 

To investigate spatial differences in the weak silicate emis- 
sion we considered the shape of the feature in the Spitzer, MIDI, 
and difference spectrum in a consistent manner. We approxi- 
mated a continuum with a straight line intersecting the spectra at 
8.3 and 13.2 fj.m, subtracted this from the spectra and then nor- 
malized the spectra with the flux level at 10.3//m. The bottom 
panel of Fig.|5] shows the result. The shape of the silicate fea- 
ture inside and outside a radius of ~ 2 AU is very similar, which 
implies that the silicate composition and grain size distribution 
should also be very similar. This is not surprising, since the high 
visibility indicates that the silicate emission is dominated by a 
compact inner region (r < 10 AU). 

4.7. Observational picture 

Before we describe a detailed modeling effort of the circumstel- 
lar material, we summarize the analyses of the various observa- 
tions. This already gives an insight into the spatial distribution of 
the gas and the dust. In Fig.|6]we display schematic representa- 
tions of two disk models, that will be considered in Sect.|5] In the 
second model we indicated the diagnostic scales of the various 
data sets to combine them into a consistent picture. 

The AMBER data probes the very inner parts of the disk. The 
analysis shows that the K-band emission could be explained with 
an emitting ring at ~Q.4 AU. A more physical model, which will 
be presented in Sect.|5]has the inner rim of the disk at 0.55 AU. 
A part of the K-band emission was shown to come from close to 
the star, and well within the inner dust disk (represented by the 
ring). This could be an indicator of ongoing accretion. 

The [OI] 6300 A emission line is form ed by the photo - 
dissociation of OH molecules by UV photons (Acke et al.l2005h . 
The detection of the [OI] 6300 A at large distances (from one to 
tens of AU) from the central star is thus an indication that the 
outer disk has an illuminated gas surface. 

The Spitzer data establish the presence of PAH emission 
and the VISIR spectrum pins it down to a circumstellar disk. 
The emission features of PAH molecules are caused by inter- 
nal vibrational modes, which are mainly excited by UV photons. 
Since the PAH molecules are coupled to the gas, the PAH emis- 
sion is another indicator of an illuminated gas surface. The re- 
solved VISIR spectrum sets the radial scale of this gas surface at 
-lOOAU. 

The VISIR Q-band image displays a faint extended axi- 
symmetric emission component (~10% of the total flux), that 
stretches out to large radii (r ~150 AU). 
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The MIDI correlated flux spectrum shows that the PAH fea- 
tures originate from radii much larger than ~2 AU and that the 
silicate composition is quite similar in the inner and outer disk. 

5. Modeling 

In this section we present a model for the disk around HD 9588 1 
which quite accurately reproduces the observations described 
above. To obtain this model we^ use the Monte Carlo radia- 
tive transfer code MCMax bv lMin et al.l (12009 ). This code can 
compute a self-consistent disk structure and a full range of ob- 
servables. It has a build-in option that models the full PAH ex- 
citation, using the temperature distribution approximation (see 
e.g. Guhathakurta & Draine 1989; Siebenmorgen et al. 1992) in- 
cluding multi-photon events for the excitation. We use MCMax 
here to compute the temperature structure, the vertical den- 
sity structure and the resulting SED, the Spitzer spectrum, the 
AMBER visibilities, the MIDI correlated flux, the VISIR im- 
ages, and the FWHM as function of wavelength. The steps to 
come to the final model presented here were the following. 

5. 1 . Initial constraints 

First we fixed the composition, and the size and shape distri- 
bution of the silicate component of the dust t o be equal to that 
obtained from the 10 micron silicate feature bv lvan Boekel et all 
(l2005h . For HD 95881 these are 80% large (1.5^m) pyroxene 
grains, 11% large enstatite grains, 1% small (0.1 yum) forsterite 
grains, 5% large forsterite grains, and 3% small (0.1 jum) sil- 
ica grains. For the shapes of t he grains we use the Distribution 
of H ollow Spheres (DHS, see iMin et al.L l2005h . as outlined in 
Ivan Boekel et al.. (.2005) . where we refer to for further details on 
the model for the silicate dust grains. In order to get the required 
continuum opacity needed we added amorphous carbon grains. 
To model the irregular shape of the carbon grains we also used 
DHS. F or the refractive index of carbon we adopted the data 
by Preibisch et alj (Il993h . Note that the continuum component 
is most likely not all in the form of amorphous carbon. Small 
grains of metallic iron and/or iron sulfide have extinction prop- 
erties similar to that of carbon. Also, large grains of various dust 
species could produce the observed continuum component. The 
abundance of amorphous carbon is a fitting parameter We have 
to realize that besides this observable dust component, a 'hid- 
den' component of large grains probably exists deep inside the 
disk. Information on this dust component, which might contain 
most of the dust mass, can only be obtained through millimeter 
observations which are currently not available for this source. In 
the following we will use the term 'visible' for the dust grains 
we can address with the current set of observations. 

The other free parameters all have to do with the geometry 
of the disk. As a first step we focused on the thermal dust grains, 
ignoring the PAH bands. The density distribution of the d ust disk 
was p arameterized using a radial surface density (Hughes et alJ 
120081) 

^Ocxr-^expl-J^j "j, (1) 

for Rin < r < Roui- Here /?() is the turnover point from where an 
exponential decay of the surface density sets in and p sets the 
powerlaw in the inner region. We fixed this powerlaw to p = 1, 
a commonly used value (see e.g. Dullemond et al. 2006). Using 
p - I basically implies that we assume the viscosity of the disk 
to linearly increase with radius. For a discussion on the possible 



values of p and their implication s for the origin of t urbulence and 
viscosity in the disk we refer to Isella etal T (l2009l) . The vertical 
density structure was computed from hydrostatic equilibrium. 
However, since we found that this vertical structure cannot re- 
produce the SED, we included a scaling parameter ^ by which 
the scale-height of the entire disk is increased. For the density 
distribution we thus have the inner and outer radii Rin and 7?out, 
the turnover radius /?o, the total dust mass in visible grains Mdust, 
and the vertical density scaling parameter ^', as free parameters. 

5.2. Fitting fhe SED 

In order to get the observed peak over continuum ratio of the 10 
and 20 micron silicate features observed in the Spitzer spectrum, 
we find that the mass in amorphous carbon grains has to be 25% 
of the total dust mass. Again, we stress that this component rep- 
resents all possible sources of continuum opacity, including very 
large grains. Second, we constrained the parameters describing 
the density structure. We found that the vertical height of the 
disk needs to be significantly increased compared to hydrostatic 
equilibrium in or der to obtain the large near IR excess (see also 
lAcke et aDl2009l for a discussion on the inner rim height). The 
height of the disk was scaled with a factor - 2.75. The ex- 
ponential decay of the surface density sets in at Rq = 2.5 AU. 
To put the density structure in some perspective, at -7.5 AU our 
prescription produces the same density as does a 1/r^ density 
law, but beyond 10 AU the surface density becomes negligible. 
The inner radius of the disk is at 0.55 AU, while the outer ra- 
dius has no influence on the observationally constrained part of 
the SED as long as it is beyond ~10AU. We find a total mass 
in visible dust grains of 10"^ Mq. Note that the total dust mass 
is probably much higher, because the mass in large dust grains 
is not constrained. To this dust distribution we added PAHs in 
an abundance needed to explain the features seen in the Spitzer 
spe ctrum. For the opac ity of the PAHs we use those computed 
bv Draine & TH (1200 lb for molecules consisting of 100 carbon 
atoms. 

5.3. Spatial distribution of tiie gas 

We consider two possibilities for the spatial distribution of the 
gas, as traced by the PAHs (see Fig.|6]l. The first is to assume 
that the gas and the dust have the same spatial distribution in 
the disk. The second is to assume that the gas does not have 
the exponential decay of the surface density for radii larger than 
2.5 AU, but that this disappearance of the dust at these radii is 
caused by grain growth and settling, which do not affect the gas. 

The first model, which has the gas and dust in the same 
spatial distribution, results in a fairly large abundance of PAH 
molecules needed to explain the strength of the features seen. 
The PAHs are in this case fairly well shielded from the stellar 
radiation by the dust grains, and thus a large amount is needed. 
Furthermore, in order to explain the absence of PAH emission 
in the MIDI correlated flux, we find that we have to remove 
the PAHs from the inner 2 AU. The total mass in PAHs in this 
case is 3-10"^ M©. This model does not reproduce the increase 
in FWHM at the wavelength of the PAH features that we found 
in the VISIR spectroscopic data (see Sect. 14. 5b . 

The second model, which assumes the gas is distributed in 
a more extended flared disk (see Fig.|6]l resulted in a very good 
overall fit of all the observables presented above. In this model 
we put the PAHs in a disk with a similar surface density as the 
dust grains but with Rq - oo, i.e. the surface density remains a 
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Fig. 6. Schematic depiction of the two considered disk models. First we consider a spatial distribution that is the same for both the 
dust and the gas. Second we consider a model that assumes that the surface density of the gas decays much slower with distance 
from the star. In this second model we indicated the different regions of the disk that are probed by the observations discussed in 
this paper 



Table 2. The parameters for our final model for the circumstel- 
lar disk of HD 95881. The carbon abundance was found to be 
25% of the visible dust. The visible dust grains and the PAH 
molecules have a different spatial distribution. 



Parameter 


Dust 


PAHs 


Inner disk radius (R,„) 


0.55 AU 


0.55 AU 


Outer disk radius (Rom) 


200 AU 


200 AU 


Inclination angle (i) 


55° 


55° 


Position angle (major axis E of N) 


103° 


103° 


Vertical density scaling parameter C¥) 


2.75 


2.75 


Power law for the surface density (p) 


1 


1 


Turnover point (Rq) 


2.5 AU 


oo 


Mass (M) 


10-* Mo 


5- 10' Mo 



powerlaw (/? = 1) for all radii. In this way we create a flaring 
outer gas disk, which is able to catch much of the radiation from 
the star. We find that the PAH emission is dominated entirely by 
flie outer regions. If the PA Hs would be destroyed according to 
the mechanism proposed bvl Dullemond et al] ( l2007h we find that 
the inner 25 AU should be free of PAHs. However, even without 
destruction the contribution to the PAH emission from these in- 
ner regions is negligible. Thus, we cannot confirm whether or not 
PAH destruction takes place in this disk. Extrapolating the pow- 
erlaw surface density distribution of the PAHs to the inner edge 
(at 0.55 AU) we find that at Rin the PAH abundance is 0.25% of 
the visible dust mass. The total PAH mass we find is 5-10-^ M©, 
which is large compared to the mass in visible dust grains, be- 
cause the PAH disk is so much larger. Including PAH destruction 
in the inner disk would only lower the total PAH mass by 12%. 
Note, that the total dust mass, including the well shielded mid- 
plane grains, can be much higher than the mass in visible grains 
and the PAH mass. This could be confirmed by future millimeter 
observations. 



5.4. Final model 

The second model, that assumes that the PAHs do not have an 
exponential decay with radius, reproduces the available spatial 
information much better This second model was fine-tuned to 
best reproduce all available observables. The position angle and 



inclination of the disk were constrained using the interferometric 
data. We find that the AMBER interferometric observations put 
important constraints on the parameters for the inner regions of 
the disk, lifting some of the degeneracies present when these 
data are not considered. The parameters of the final model are 
summarized in table|2] 



6. Discussion 

6.1. Comparison of the model with the observations 

Our disk model accurately reproduces all the observables. The 
fits to the SED, the Spitzer spectrum, the MIDI correlated flux, 
the AMBER visibilities, and the VISIR spectroscopic FWHM 
are presented in Fig.|7] 

Spitzer. The general slope and most features in the Spitzer 
spectrum are reproduced. However, the observed 6.2, 7.9 and 
8.6 yum PAH features are stronger than in our model and the pre- 
dicted PAH features around 20jum are not seen in the Spitzer 
spectrum. These differences are related to the PAH chem istry, 
which is a much debated subject (see e.g. i Draine & Lill2007l) . 

MIDI. The N-band correlated flux obtained by MIDI has an 
error on the absolute calibration of approximately 10%, similar 
to the diff'erence with the model output. Thus our modeling of the 
spatial distribution of the visible dust grains is consistent with 
the MIDI result. 

AMBER. The simple flat ring+point-source model of 
Sect.|4T| actually reproduced the inclination, position angle and 
visibilities of the disk very well, although the inner rim radius 
(0.37 AU) is significantly below the 0.55 AU of our final model. 
Our more physical final model gives a slightly poorer fit to the 
visibilities. This indicates that the structure of the inner disk is 
more complicated than assumed. The exact structure is currently 
much debated (see Sect. 16. 21 ). 

VISIR spectrum. The model FWHM was obtained by making 
a Gaussian fit to the spatial profile after convolving the model 
with a Gaussian of the same width as the PSF of the VISIR 
11.2yum setting. The PSF of the 8.5 jum setting was scaled to 
the model. The matching continuum levels are thus a result of 
our method, but the increase in the FWHM at the PAH wave- 
lengths are a confirmation of the correct modeling of the spatial 
distribution of the PAH emission. 
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Fig. 7. The best model fit to the observables of HD 9588 1 . The upper left panel shows the fit to the spectral energy distribution with 
an inset for the Spitzer spectrum. The blue line is the observed spectrum, the black line the model spectrum and the points with 
error bars are photometric measurements obtained from the literature. The upper right panel shows the correlated flux as obtained 
by MIDI. The blue line is the observed correlated flux and the black line the model. The bottom left panel shows the visibilities 
obtained by AMBER. The colored lines give the model results (color coding the same as Fig.[T]). The bottom right panel shows in 
blue the VISIR FWHM as a function of wavelength and in black the model. 



VISIR image. A model of the Q2 image was obtained by tak- 
ing the output of our model at 19.0 fim, just next to the artificial 
PAH feature. We convolved this model image with the VISIR 
PSF and we applied a PSF subtraction in analogy to the image 
analysis of Sect. l4.4l The observed resolved emission compo- 
nent is not reproduced by the model. A likely explanation for 
this is the photoluminescence of particles with sizes between 
dust grains and PAH molecules, the so c alled very small g rains, 
which are not included in our model (see lFlagev et al.ll2006i) . 



6.2. The distribution of gas and dust 

Our modeling reveals that the spatial distribution of the gas and 
the dust in the disk around HD 95881 is different. The dust in 
the upper layers of the outer disk is heavily depleted, while 
the gas still has a large scale height. To reproduce the near IR 
flux we needed to assume a scale height, which is much larger 
than would be obtained from vertical hydrostatic equilibrium 
(*? = 2.75). We do not have an explanation for this, but we 
note that the value of is sensitive to the assumptions on the 
structure and composition of the material in the inner disk (see 



e.g. Isella et al., 2008; Tannirkul am et al ],L2008aV Furthermore, 
our value is well withi n the range of scaling parameters (1.0-3.0) 
needed bv lAcke et al.l (|2009) to explain the near IR SED of about 
half of their sample of ~30 Herbig Ae/Be stars. 

The spatial distribution of the visible grain component was 
well constrained by the SED. The large grains could not be con- 
strained directly, but are probably abundantly present in a settled 
outer disk. Observations of the millimeter flux will help to con- 
strain the mass in this component as well as the size distribution. 
We can however make an estimate of the total dust mass, based 
on the modeled total PAH mass (5 10 "^M©) and the modeled 
PAH abundance at the inner edge (0.25%). Assuming that at the 
inner rim all the dust is visible and assuming that the PAH to 
dust ratio is homogeneous throughout the disk we find that the 
total dust mass adds up to 2-10"*Mq. This means that most of 
the dust mass resides in large grains (200: 1). Using the canonical 
gas to dust ratio of 100 the total disk mass becomes 2-10"'* Mq. 
Note, that if the outer disk has a lower value of "F, a larger gas 
mass could be required to reproduce the observations. 

The results presented above naturally lead to the picture of 
a disk in which the dust grains in the outer disk are coagulated 
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and settled towards the midplane, while the gas is still available 
to keep the PAH molecules in the higher atmosphere of the disk. 
As alr eady noted in the theoretical study by Dullemond et al. 
(I2007h . growth and settling of the dust grains leads to a natu- 
ral increase of the relative strength of the PAH signature, as is 
observed in this disk and confirmed by our modeling effort. 

6.3. Context 

In general the study of iDullemond et al.l (l2007h showed that 
the natural outcome of a group I source after grain-growth and 
sedimentation of the dust is a group II source, that maintains 
the flaring structure for the gas. However, observational studies 
have sho wn that most group II sources l ack a flari ng gas dis- 
tribution. iMeeus et al.1 (l200ll) . lAcke et al.1 (|2004 and lAcke etal] 
(1201 Oh showed that group I sources display significantly more 
PAH emission. lAcke et al.1 dlool showed that group I sources 
have in general stronger [OI] emission. Apparently the gas 
of most group II sources has either dramatically decreased its 
scale height because of the lack of heating or the gas has been 
dispersed from the disk. How disks can lo se their gas is cur- 
rently being debated (see lHillenbrandll2008h . but photoevapora- 
tion seems to be the most likely mechanism. 

On the other hand there is a fair number of group II 
sources that do show indications of a flaring gas distribu- 
tion. Some group II sources displ ay PAH emission in their 
1 yum spectr a, for example: HKOri jvan Boekel et al I2OOI and 
HD 142666 ( Verhoefi'et al.ll2010 !). Some group II sources dis- 
play the [O I] 630 A Une in emission, for example: HD 98922 
(lAcke et al]|2005h . And some sho w bofli PAH a nd [0 11 emis- 
sion, for example: HD 101412 (" Fedele et al] l2008l) . These 
sources all seem to be in a transitional phase from a gas rich 
flaring disk to a gas poor self-shadowed disk. Roughly half of 
the known group II sources are in this phase, which means we 
can infer that half of the life time of the disk of a group II source 
is spent on the dispersal of the gas. An estimate for this time 
scale is t he photoevaporation time s cale, which is on the order of 
~ 1 0^ yr (iGorti & Hollenbachll2009h 

7. Conclusions 

A comprehensive study was performed to map the distribution of 
the gas and dust in the protoplanetary disk around HD 95881. In 
the right panel of Fig.|6]we displayed a schematic representation 
of the disk, which puts all results in perspective. The AMBER K- 
band interferometry showed that there is an extended hot inner 
region with emission coming from well within the sublimation 
radius. The detection of the [O I] 6300 A indicated that the disk 
has a flaring gas surface at large distances (from one to tens of 
AU) from the star. The finding of PAH features in the Spitzer 
and VISIR spectra confirmed the presence of an illuminated gas 
surface. The resolved VISIR spectrum traced this surface up to 
radii of ~ 1 00 AU. 

W e used the radiative transfer code MCMax (iMin et alJ 
l2009h to create a model of the disk's density and temperature 
structure. Our model satisfactorily reproduced all of our obser- 
vations. The main conclusions that followed from our model are 
that the inner disk contains most of the visible grains and has a 
puffed up inner rim, the dust grains in the outer disk have co- 
agulated and settled towards the midplane, while the gas and 
PAH mixture maintain a flaring geom etry. Theory predic ted the 
existence of these type of disks (Dul lemond et al.ll2007i) . while 
observational trends showed that most of the sources with self- 



shadowed dust distributions have dispersed the bulk of their gas. 
In this light HD 95881 is a special source: it is in the transi- 
tion phase from a gas rich flaring dust disk to a gas poor self- 
shadowed disk. 
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